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ABSTRACT 

Aims. Our aim is to perform on galaxy mock catalogues the same colour-density analysis earned out by Cucciati et al. (2006) on a 5/r'Mpc scale 
using the VIMOS-VLT Deep Survey (VVDS), and to compare the results from mocks with observed data. This allows us to test galaxy evolution 
in mocks, and to understand the relation between the studied environment and the underlying dark matter distribution. 

Methods. We used galaxy mock catalogues, with the same flux limits as the VVDS-Deep (Iab ^ 24) survey (Cmocks), constructed using the 
semi-analytic galaxy catalogues by De Lucia & Blaizot (2007) applied to the Millennium Simulation. From each mock, we extracted a sub-sample 
of galaxies mimicking the VVDS observational strategy (Omocks). We then computed the B-band Luminosity Function LF and the colour-density 
relation in mocks using the same methods employed for the VVDS data. 

Results. We find that the B-band LF in mocks roughly agrees with the observed LF, but at 0.2 < z < 0.8 the faint-end slope of the model LF is 
steeper than the observed one. Computing the LF for early and late type galaxies separately, we show that mocks have an excess of faint early 
type galaxies and of bright late type galaxies with respect to data. We find that the colour-density relation in Omocks is in excellent agreement 
with the one in Cmocks. This suggests that the VVDS observational strategy does not introduce any severe bias to the observed colour-density 
relation. At z ~ 0.7, the colour-density relation in mocks agrees qualitatively with observations, with red galaxies residing preferentially in high 
densities. However, the strength of the colour-density relation in mocks does not vary within 0.2 < z < 1.5, while the observed relation flattens 
with increasing redshift and possibly inverts at z ~ 1.3. We argue that the lack of evolution in the colour-density relation in mocks cannot be due 
only to inaccurate prescriptions for the evolution of satellite galaxies, but indicates that also the treatment of the central galaxies has to be revised. 
Conclusions. The reversal of the colour-density relation can be explained by wet mergers between young galaxies, producing a starburst event. 
This should be seen on group scales, where mergers are frequent, with possibly some residual trend on larger scales. This residual is found in 
observations at z = 1.5 on a scale of ~5/r'Mpc, but not in the mocks, suggesting that the treatment of physical processes influencing both satellites 
and central galaxies in models should be revised. A detailed analysis would be desirable also on small scales, which requires flux limits fainter 
than those of the VVDS data. 
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1. Introduction 

Several physical mechanisms are expected to influence the 
properties of gal axies in over-dens e regions: ram pressure 
stripping of gas (|Gunn & Gottl Il972h . gal axy-galaxy mergin g 
dToomre & ToomreT 1972b. strangulation (lLarson et al.l | 1980l) . 
and harassment (iFarouki & Shapirolll98ll iMoore et al.lll996l) . 
Each mechanism has specific environmental dependencies and 
timescales, but their relative role in regulating galaxy formation 
and determining the observed trends remains unclear. In addi- 
tion, it is not yet clear to what extent internal processes such as 
feedback from supernovae and central black holes contribute to 
the observed environmental dependence of the galaxy structural 
parameters. Finally, it is known that in a Gaussian random field 
there is a statistical correlation between mass fluctuations on 
different scales, with most massive h alos preferen t ially residing 
withi n large scale over-densities (see lKaiserin~987t iMo & White! 
119961) . What is less clear is the role of initial cosmological condi- 
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tions in m odulating the observed density dependence of galaxy 
properties (lAbbas & Shethl2005h . 

Semi-analytic models (SAMs) of galaxy formation coupled 
with dark matter (DM) simulations provide a useful tool to ad- 
dress these issues. In SAMs, the individual physical processes 
taking place during galaxy evolution are expressed through sim- 
ple equations that are motivated by observational and/or theo- 
retical studies. These equations parametrise the dependency of 
such processes on physical properties of galaxies and/or those 
of the dark matter haloes in which they reside. Then, the com- 
parison between galaxy properties in SAMs and in the observed 
data gives important feedback on the validity of the prescriptions 
used in the models and on the nature of the observed correlation 
between galaxy properties and environment. Models are usually 
tuned to reproduce some (sub)set of observational data in the lo- 
cal Universe, most notably the observed local galaxy luminosity 
function (LF) and/or mass function. Tracing back in time galaxy 
evolutionary paths, pushing these studies at high redshifts, be- 
comes thus a very powerful tool in terms of discriminating be- 
tween different models. 



1 



O. Cucciati et al,; The VVDS colour-density relation at z < 1.5: model vs observations. 



The recent completion of large and deep high-redshift sur- 
veys makes a detailed comparison betwe en model pred i ctions 
and observational dat a possible (see e.g. IStringer et al.l 120091 : 
Ide la Torre et al.ll2.01 ll) . Following galaxy evolution over a large 
redshift range can help addressing the open issue of the apparent 
contradiction between the hierarchical growth of DM structures 
and the 'downsizing' scenario of luminous matter. A number of 
different observational tests support a 'hierarchical' scenario for 
structure formation in which DM halos form first in the high- 
est density peaks of the primordial density field, and then grow 
hierarchically through subsequent mergers (e.g. lPeebleslll98Cib . 
Observations have shown that galaxy evolution does not proceed 
in a similar 'bottom-up' fashion, at least for their star formation 
(SF) histories: with in creasing cosmic time, SF moves toward s 
less massive galaxie s (ICowie et al. 1996t Gavazzi et a 
As discussed in e.g. iDe Lucia et al. (l2006hT these findings are 
not necessarily evidence of 'anti-hierarchical' growth of lumi- 
nous matter as the 'formation' of the gala xy stellar population 
does not coincide with its assembly (see iFontanot et al.l 120091 
for a more detailed discussion about different manifestations of 
downsizing and comparison with observational data). 

Simulations can also help clarifying the relation between dif- 
ferent definitions that are commonly adopted for the 'environ- 
ment' . Indeed, different quantities have been used in the litera- 
ture to characterise the local and/or global environment. The use 
of different environmental definitions makes it very difficult to 
compare results from different surveys, and at different cosmic 
epochs. Using simulated galaxy catalogues allow us to have a 
common reference for the environment, such as the underlying 
DM distribution. 

In this paper, we use mock catalogues constructed from 
semi-analytic model s applied to the Millennium Simulation^ 
dSpringel et al.l 120051) to carry out a detailed comparison with 
the observational results presented in ICucciati et al.l (1200 6l C06 
her eafter), based on the VIMOS-VLT Deep Survey (VVDS, 
see iLe Fevre et al. 1 12001 . A detailed comparison between ob- 
servations and model data for galaxy number counts, redshift 
and colour distribution , and galaxy clustering was presented in 
Ide la Torre et al.l (1201 ll Paper I hereafter). Here we focus on the 
observed colour-density relation and its evolution. 

The aims of our study are: (i) test the robustness of observa- 
tional results in C06 versus possible biases due to the observa- 
tional strategy adopted (e.g., does the VVDS sampling rate alter 
the strength of the colour-density relation?); (ii) test galaxy evo- 
lution in mocks (do we see the same environmental effects on 
galaxy properties in mocks and observed data?); (iii) understand 
what 'environment' is being studied (e.g., which is the relation 
between the density traced by galaxies and the underlying DM 
distribution?). 

This paper is organised as follows. Section 2 gives a sum- 
mary of the VVDS data used in C06 and describes the mock 
galaxy catalogues used for the model comparison. In Section 3, 
we compare the rest-frame B-band LF of the VVDS data with 
that obtained using mock catalogues. Section 4 describes the lo- 
cal density and colour distributions in mocks, and compares the 
colour-density relation in C06 with the one found in mocks. In 
Section 5 we take advantage of the available mocks to analyse 
how the local density computed on a 5/i _1 Mpc scale compares 
with the total halo mass in which galaxies reside. In Section 6 
we discuss our results, and we summarise them in Section 7. 

Throughout this paper, we use the AB flux normalisation for 
both observed data and mocks. When we refer to observed data, 



we adopt the concordance cosmology (Q,„, Qa, h) = (0.3, 0.7, 
0.7). 



2. Data and mock catalogues 

2.1. The VVDS Deep sample 

The VVDS is a large spectroscopic survey, with the primary aim 
of studying galaxy evolution and large scale structure formation 
over the redshif t range < z < 5. Det ails on the survey strategy 
are described in lLe Fevre et al.l (I2005I) . VVDS is complemented 
by ancil lary deep photometric data: BVRI from the C FHT-12K 
camera dMcCracken et al.ll2003ULe Fevre et al.ll2004l) JK f rom 
the NTT telescope dlovino et alj|2005t iTemporin et al.1l2008l) . U 
from the MPI telescope dRadovich et al.ll2004l) . u*. s'. r\ V, z'- 
band data from the CFHT Legacy Survey and JHK$ from the 
WIPvDS sur vey dBielbv et al.1 ifoTll) with the CFHT-WIRCAM 
camera (see lCucciati et al.ll2012 for a detailed description). 

This paper is based on the colour-density relation stud- 
ied in C06, over the VVDS-0226-04 Deep field (from now on 
"VVDS-02h field"). We refer the reader to that paper for a 
detailed description of the data. Briefly, the VVDS-02h data 
set is a purely flux limited spectroscopic sample, with 17.5 < 
Iab ^ 24.0. In this range of magnitudes, the parent photomet- 
ric catalogu e is complete and free fr om surface brightness selec- 
tion effects dMcCracken et al.l2003l) . Spectroscopic observations 
were carried out at the ESO-VLT with the Visible Multi-Object 
Spectrograph (VIMOS) using the LRRed grism. The rms accu- 
racy o f the redshift measurements is ~ 275 km/s (ILe Fevre et al.l 
120051) . The VVDS-02h field covers a total area of 0.7 x 0.7 deg 2 , 
targeted by 1, 2 or 4 spectrograph passes, and it probes a co- 
moving volume (up to z=1.5) of nearly 1.5 x 10 6 /i~ 3 Mpc 3 in 
a standard ACDM cosmology. The covered field has transver- 
sal dimensions ~ 37x37 /i~'Mpc at z=1.5. Averaging over the 
area observed, spectra have been obtained for a total of 22.8% 
of the photometric sources, and ~ 80% of these targeted objects 
yield a relia ble redshift, resul ting in an overall sampling rate of 
~ 18% (see lllbert et al.ll2005l) . The final galaxy sample used in 
C06 contains 6582 galaxies with reliable redshifts. 



http://www.mpa-garching.mpg.de/galform/virgo/millennium/ 



2.2. The mock catalogues 

Mock galaxy catalogues were obtained by applying the 
semi-analytical model o f galaxy evolution described in 
IDe Lucia & Blaizol d2007l) to the dark mat ter halo merging tree s 
derived from the Millennium Simulation dSpringel et al.l 120051) . 
This contains N = 2160 3 particles of mass 8.6 x 10 8 /i _1 M o 
within a comoving box of size 500 /j _1 Mpc on a side. The 
adopted cosmological model is a ACDM model with Q.„, = 0.25, 
Q b = 0.045, h = 0.73, fi A = 0.75, n = 1 and cr & = 0.9. 

The semi-analyt i cal mo del used here is fully described in 
IDe Lucia & Blaizol ( 2007[). It builds on results from previ- 
ous works dKauffmann & Haehneli 120001: ISpringel et~aDl200lb 
IDe Lucia et alJl2004tlCroton et al.ll2006l) to describe the relevant 
physical processes, e.g. star formation, gas accretion and cool- 
ing, formation of super-massive black holes and galaxy mergers, 
feedback (also 'radio mode' feedback). The SAM model used in 
this study has been tested against many observational data and 
and it has been shown to provide a relatively good agreement 
with observational data both for the local Universe and at higher 
redshift. It is, however, not without problems. In particular, the 
fraction of low-mass red galaxies is too high, and the clustering 
signal of red galaxies is overpredicted. The galaxy mass func- 
tion has a too high normalisation in the regime of low and in- 
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ter mediate stellar masses, a t any redshif t s. We r efer the reader 
to IWeinmann et all (120061) . IWang et al.l (120081) . iFontanot et aD 
d2009t) and Paper I for more details. It is worth noticing that 
these problems are not specific of the particular model used here 
but appear to be common to most recently published SAMs. We 
will further discuss some of these discrepancies in the following 
sections. 

For our analysis, it is important to summarise how abso- 
lute magnitudes are com puted in the SAM adopted here. The 
iBruzual & Charlotl (120031 ) model is used to generate lookup ta- 
bles for the luminosity of a single burst of fixed mass, as a func- 
tion of the age of the stellar population (f) and its stellar metal - 
licity (Z). At each star formation episode, the model interpolates 
between these tables, using a linear interpolation in t and log(Z), 
to calculate the contribution to the luminosity of model galax- 
ies at the time of the observations. Stars are assumed to form 
with the metallicity of the cold gas component, and an instanta- 
neous recycling approximation is adopted. The model adopts a 
IChabrierl (120031) IMF with a lower and upper mass cut-off of 0. 1 
and 100 M Q , respectively. Magnitudes are e xtincted (internal ex- 
tinctio n) using the model that is detailed in lDe Lucia & Blaizotl 
J2007h . 

For our study, we randomly selected 20 1 x 1 deg 2 
Millenni um light cones, fro m those constructed using the code 
MoMaF dBlaizot et a l. 2005) as described in Paper I. From these 
cones, we extracted two sets of mock catalogues. First, we ex- 
tracted lxl deg 2 mocks with the same flux limits as the VVDS- 
02h sample (17.5 < Iab ^ 24), and added a posteriori the same 
redshift measurement error of the VVDS sample. These cata- 
logues, having 100% sampling rate, will be called from now on 
Cmocks. Then, from the Cmocks we extracted mock catalogues 
mimicking our VVDS-02h sample, i.e. adding geometrical ef- 
fects and uneven sampling rate. We call these mocks Omocks. 
We refer the reader to Paper I for details. 

As we will show later (Sec. 12.2. 11 1. Millennium light cones 
have on average higher /-band number counts (per surface area 
and over the range 17.5 < I AB < 24) than VVDS. In this study 
we complement the 20 random light cones with two additional 
cones with the lowest number counts (cones 036d and 045b), 
among those presented in Paper I. VVDS number counts are in 
between the number counts of these two cones. We also added 
the cone with the highest number counts (cone 022b). 
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Fig. 1. Galaxy number counts per unit magnitude and per square 
degree. Red triangles are for VVDS measurements, with er- 
ror bars representing the Poissonian uncertainty. The red line 
shows the median number counts of the 23 mocks, and the grey 
area shows the 16''' to 84' /! percentile range of the mock num- 
ber counts distribution. The VVDS number c ounts are in good 
agreem ent with number counts obtained by iMcCracken et al.l 
(120031) for the CFH 1 2K- VIRMOS deep field. 
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2.2.1 . Mock catalogues number counts and mean 
inter-galaxy separation 

Since we will compute the galaxy 3D local environment in 
mocks, and compare it to the one found in the VVDS, it is impor- 
tant to know possible similarities and/or discrepancies between 
the 3D galaxy distribution in the mocks and that in the data. 

In Paper I, we have shown that the light-cones used in our 
study give an average redshift distribution n(z) and galaxy clus- 
tering that are not in perfect agreement with the VVDS-02h mea- 
surements. In particular, we showed that the mock n{z) agrees 
with the observed one in the range 0.5 < z < 1.8, but it overesti- 
mates the observed number of galaxies at z < 0.5, and underes- 
timates it at z > 1 .8. In addition, we have shown that the Munich 
semi-analytical model overestimates the VVDS clustering, and 
that this is mainly due to an excess of the clustering signal for 
model red galaxies. 

Fig-fflshows the number counts per square degree of VVDS- 
02h galaxies as a function of the observed Iab magnitude (that 
is the VVDS selection magnitude). Number counts have been 
computed from the parent photometric catalogue of ~40,000 



Fig. 2. Mean inter-galaxy separation as a function of redshift, in 
the four redshift bins used in C06 (0.25-0.6, 0.6-0.9, 0.9-1 .2, 1 .2- 
1.5). Red squares are the mean among the 23 Omocks used in this 
study, and the grey shaded area shows the 1 - a scatter. Blue cir- 
cles are for VVDS data. The thick solid line corresponds to the 
mock with mean inter-galaxy separation closest to the VVDS 
measurements at all redshift bins considered. The dashed line 
corresponds to the mock with lowest number counts (in the to- 
tal redshift range 0.25-1.5), and the dotted line to the cone with 
highest number counts. 



sources, with stars rem oved using the method described in 
IMcCracken et al.l d2003l) . The median number counts of the 23 
mocks used in this study are over-plotted as a solid line. Fig. Q] 
shows that the mock /-band number counts are higher than ob- 
servational data at Iab > 22, but lower than the observed number 
counts at the brightest magnitudes (Iab <20). Since faint galax- 



3 



O. Cucciati et al,; The VVDS colour-density relation at z < 1.5: model vs observations. 



ies are much more numerous than bright ones, mocks contain on 
average ~ 10% more galaxies than observations^. 

Fig.|2]shows the mean inter-galaxy separation in the VVDS- 
02h field, as a function of redshift. We overplot the same quantity 
computed in the Omocks. In the lowest redshift bin, all mocks 
considered have mean separation lower than measured in ob- 
served data. The average mean separation in mocks is larger than 
the VVDS one in the range 0.6 < z < 1 .2, and lower than the 
measured VVDS value for z > 1.2. At z > 0.6, there is at least 
one mock with mean separation close to the VVDS one, but this 
is not the same mock at all redshifts. In Fig. [2] we also show 
the mean inter-galaxy separation for the mock with the lowest 
number counts (dashed line). Its mean inter-galaxy separation is 
larger or very similar to the VVDS one (at least up to z = 1.2), 
as one would expect. The opposite is true for the mock with the 
highest number counts (dotted line), up to z — 1 .2. Among those 
used in this study, there is no mock with number counts close to 
the observed ones at all redshift bins. The black solid line in the 
figure shows the mock that has inter-galaxy separation on aver- 
age the most similar to VVDS at all redshift bins. Its counts are 
very close to the average among all the mocks at each redshift 
(red squares). It does not correspond to the mock with the lowest 
total number counts. 

The results on the n(z), the number counts as a function of 
selection magnitude, and the mean inter-galaxy separation are 
all related, showing that Omocks contain more galaxies than the 
VVDS-02h field. 



2.2.2. The absolute magnitudes 

In C06, red and blue galaxies are defined using the rest- 
frame colour u* — g' . The u* and g' (CFHT-LS filters) rest- 
frame magnitudes are not available on the Millennium database. 
Among the rest-frame ab solute magnitude available for the 
iDe Lucia & Blaizotl (120071) . the closest to u* and g' are the ones 
in the B and V Johnson filters. To be consistent with the anal- 
ysis discussed in C06, we computed for our mocks the u* and 
g' absolute magnitudes using the same method employed for 
the VVDS data (see below). We also computed for our mocks 
the absolute magnitudes in the B and V Johnson filters, so as to 
compare them with the intrinsic ones available from the simula- 
tions. In this way, we can verify that the rest-frame magnitudes 
computation does not affect the intrinsic luminosity distribution. 

For the VVDS data set, we computed the absolute magni- 
tudes using the code ALF (Algorithm for Luminosity Function, 
lllbert et al.ll2005l) . based on a SED fitting technique. To reduce 
the dependency on templates, we derive the rest frame absolute 
magnitude in each band using the apparent magnitude from the 
closest observed band, shifted at the redshift of the given galaxy. 
In this way, we minimise the applied K-correction, that depends 
on the assumed template. The set of observed magnitudes used 
to derive the absolute magnitudes in C06 included BVRI bands 
from the CFHT-12K camera, and u*, g', r', i', z'-band data from 
the CFHT Legacy Survey. We have these observed magnitudes 
in our mocks, and we used them to compute the intrinsic lumi- 
nosities in the mocks, with the code ALF. We adopted the same 



2 Fig. 3 of Paper I shows that /-band counts in the mocks are in very 
good agreement with those from the VVDS. The discrepancy with re- 
sults presented here is due to an error in the plotting routine used for 
that Figure. This however does not alter the conclusions of that section 
of Paper I, which are based on the /'-band counts (Fig. 4 in paper I): the 
/'-band counts in the mocks are higher than in the VVDS, consistently 
with our /-band results. 



cosmology used for the VVDS data, i.e. Q„, = 0.3, Qa = 0.7, 
h = 0.7. 

We compared the B and V absolute magnitudes (not in- 
cluding internal dust extinction) available on the Millennium 
database with the corresponding quantities computed with ALF. 
We found that the two sets of magnitudes are consistent within 
~ 5%. The small difference between the two sets of magnitudes 
is slightly dependent on the galaxy luminosity, and does not de- 
pend on redshift within the redshift range explored. As B and 
V filters are very close to u* and g', we are confident that also 
the computation of u* and g' rest-frame magnitudes using the 
method described above does not introduce any spurious effect. 



3. The B-band Luminosity Function 

In this section, we compare the V VDS B-band (Joh nson) 
Luminosity Func tion presented in lllbert et ail J2005) and 
IZucca etalld2006|) with the corre sponding quantity derived from 
the mocks. In lZucca et al.l d2006l) . we used only UBVR$ appar- 
ent magnitudes as input for the SED fitting, because the CFHT- 
LS photometry was then not yet available. For consistency, we 
computed a second set of absolute magnitudes for the mocks, 
using only these five observed bands. This second set of abso- 
lute magnitudes has been used only for the LF presented in this 
section. For this analysis, we have used the Cmocks catalogues, 
and compared model results with the observational estimate cor- 
rected for incompleteness. 

The LF was computed using the code ALF (see Sect. 12.2.21 . 
which i mplemen t s sev eral es timators: the non- parametric 
l/y„,„, (ISchmidll 1 19681) . C + dLvnden-Belll 1 19711). S WML 
(lEfstafhiou et al.ll 19881) . and the parametric S TY dSandage et akl 
119791). We used the S TY assuming a single Schechter function 
( Schechterll 19761) parametrised in terms of a characteristic lumi- 
nosity (L*), a faint-end slope (a), and a normalisation (density) 
parameter (<p*). For a mo re detailed descrip tion of the tool and 
the estimators, we ref er to lllbert et al.l(l2005l) . 

lllbert et all d2004l) have shown that the LF measurement can 
be biased, mainly at the faint end, when the band used is far from 
the rest frame band in which galaxies are selected. This is due 
to the fact that, because of the K-correction, galaxies of different 
type are visible in different absolute magnitude ranges at a given 
redshift, even when applying the same flux limits. Moreover, in a 
flux-limited survey, this limit varies with redshift. When comput- 
ing the VVDS LF, we avoided this bias by using in each redshift 
range only galaxies within the absolute magnitude range where 
all the SEDs are observable. We have computed the LF in mocks 
using the same method. 

The upper left panel of Fig. [3] shows the luminosity func- 
tion for one of the Cmocks catalogues (cone 001b), in differ- 
ent redshift bins, obtained with C + and STY methods. The lu- 
minosity functions derived with the other two methods (1 /V max 
and S WML) are consistent with those shown in the figure. The 
dashed red line and red filled circles in each panel show the cor- 
responding results for the VVDS data. The vertical dashed line 
represents the faint absolute limit considered in the STY esti- 
mate. 

The upper right panel of Fig. [3]shows the confidence ellipses 
of the a and M* parameters obtained in Cmock 001b (black el- 
lipse), and also in Cmock 022b (magenta ellipse) and Cmock 
036d (cyan ellipse), which are the mocks with the lowest and 
highest number counts (see Sect. 12.2b . The red thicker ellipse 



3 (/-band from the MPI telescope, and BVRI from the CFHT-12K 
camera as described in Sec. 12. II 



4 



O. Cucciati et al,; The VVDS colour-density relation at z < 1.5: model vs observations. 



refers to the VVDS parameters, and the red thick lines in the 
first and last redshift bins show the uncertainties on one of the 
parameters derived by keeping the other one fixed. This figure 
shows that shape parameters a and M* measured for the three 
mocks considered are consistent within the uncertainties. Some 
differences can be found in the normalisation parameters <p*, re- 
flecting the different number counts in the mocks. 

Overall, there is a reasonable agreement between the LF in 
the VVDS and in mocks, but with some non negligible discrep- 
ancies. In particular, we find that M* is almost always brighter 
and a almost always flatter in the VVDS than in the mocks. In 
addition, there are significant differences in the normalisation of 
the model and data LFs, reflecting the differences in the num- 
ber counts discussed in Sect. 12.2.11 In order to better explore 
these differences and to understand if they are induced by a spe- 
cific class of objects, we derived the LF for galaxies of different 
types. 

IZucca et"al] d2006l) split the global galaxy population in dif- 
ferent spectro-photometric types. For each galaxy, they found 
the best template fitting th e galaxy SEP, c hoosin g among four 
empirical templates from IColeman et al.l (Il980i CWW here- 
after) and two starburst tem plates computed using GISSEL 
dBruzual A. & Charlotlll993l) . They defined four galaxy types, 
corresponding to the four CWW templates (E/SO, early spiral, 
late spiral and irregular - type 1, 2, 3, and 4 respectively). Type 
4 galaxies include the two starburst templates. 

Here we apply the same classification s cheme to galaxies i n 
the mocks, and derive their LFs as done in IZucca et al 
The lower panels of Fig. [3] shows the LFs obtained for the two 
extreme galaxy populations (type 1 in the left bottom panel, 
and type 4 in the right bottom panel) in the mock 001b and in 
the VVDS data. Results from the mocks with largest and low- 
est number counts are similar. The figure shows a clear excess 
of type 1 galaxies in the mocks at magnitudes fainter than the 
'knee' of the LF. In contrast, bright type 1 galaxies are under- 
represented in mocks. For the latest-type galaxies (type 4), in 
mocks there is an excess of bright galaxies in the lowest red- 
shift bins considered, and a slight deficit of fainter galaxies. 
Therefore, the steeper faint-end slope and the fainter M* found 
in the global LF in mocks are due to an excess of faint type 1 
galaxies and deficit of bright type 1 galaxies, respectively. The 
effect on the global LF of the deficit of bright type 1 galaxies 
is less evident because it is compensated by the excess of bright 
type 4 galaxies. 

The excess of faint red galaxies in the model used in this 
paper is a known problem which actually appears to be shared 
by most (all) semi-ana lytic models that have been published 
in recent years (see e.g . IWang et aill2007l: iFontanot et alj|2009t 
IWeinmann et al.ll2010l and references therein). This is also con- 
sistent with what we found in Paper I. 



4. The colour-density relation 

In C06, we analysed the colour-density relation as a function of 
both redshift and galaxy luminosity. Namely, we split our sam- 
ple in four redshift bins (0.25 <z< 0.6,0.6 <z< 0.9,0.9 <z< 
1 .2 and 1.2 < z < 1 .5), and in each redshift bin we studied 
the colour-density relation for galaxies with (Mb - 5 log li) < 
-19.0, -19.5, -20.0, -20.5, -21.0. 

We stress that the density has been computed using the en- 
tire sample available, irrespectively of galaxy luminosity (see 
Sect.l4~T1i. Moreover, at 0.25 <z< 0.6, the VVDS sample does 
not contain enough galaxies with [Mb- 5 log h) < -20.5 because 



of the small probed volume. So we excluded the two brightest lu- 
minosity thresholds from the analysis in this redshift range. We 
also know that VVDS samples brighter and brighter galaxies at 
higher redshift, due to its flux limit. So we examined only galax- 
ies with {Mb — 5 log h) < -19.5 and < -20.0 at 0.9 <z < 1.2 and 
1.2 <z< 1.5, respectively. These lower luminosity limits as- 
sure that the considered sub-samples are complete for all galaxy 
types. 

4. 1 . The environment parameterisation 

We refer the reader to C06 for a detailed description of the den- 
sity computation method. Here we only give a brief summary. 

For each galaxy at a comoving position r, C06 characterised 
the environment that surrounds it by means of the dimension- 
less 3D density contrast 6(r, R), smoothed with a Gaussian filter 
of dimension R: S(r,R) = [p(r,R) -p(r)]/p(r). When smooth- 
ing, galaxies are weighted to correct for various survey observa- 
tional characteristics (sample selection function, target sampling 
rate, spectroscopic success rate, and angular sampling rate). 
Moreover, underestimates of 5 due to the presence of edges have 
been corrected by dividing the measured densities by the fraction 
of the volume of the filter contained within the survey borders. 

In C06, we calibrated the density reconstruction scheme 
using simulated m ock catalogues extracted from GallCS 
dHatton et all l2003h . The aim was to determine the redshift 
ranges and smoothing length scales R over which our environ- 
mental estimator reliably reproduced the underlying galaxy en- 
vironment, as given by a 100% sampling rate catalogue with 
Iab ^ 24. We concluded that we reliably reproduced the under- 
lying galaxy environment on scales R > 5/i~'Mpc out to z=1.5. 
In the present study, we do not use GallCS mocks because they 
do not provide the information needed for our analysis (e.g. the 
observed magnitudes in all the VVDS bands). However, for a 
sanity check, we repeated the tests carried out in C06 using the 
mocks used in the present study, and we confirmed our previous 
results on the reliability of the density reconstruction. 

We computed the density field in both the Cmocks and 
Omocks, using the same method as in the VVDS data, and the 
same flux-limited {Iab ^ 24) tracers population. As mentioned 
in Sec. l2.2l for the Omocks we have adopted the same weighting 
scheme used for the VVDS data. In particular, we have consid- 
ered the target sampling rate, the spectroscopic success rate, and 
angular sampling rate, and also accounted for survey boundaries. 
By construction, no correction is needed for the C mocks since 
they have a 100% sampling rate. Moreover, to compute the den- 
sity in the Cmocks, we started from the lxl deg 2 mocks from 
which they were extracted, so the density field depends very lit- 
tle (if any) on the correction for boundary effects. 

Fig.|4]shows the density contrast distribution for Omocks and 
VVDS data, for the different redshift bins and luminosity limits 
explored in C06. The average density distribution in Cmocks is 
very close to that in Omocks, but Cmocks have a smaller scatter 
around the mean. The density distribution in mocks has longer 
tails towards high densities than VVDS, in particular at the low- 
est and highest redshift explored. Interestingly, these longer tails 
are not consequence of the higher number counts in mocks: the 
figure shows that also the mock with lowest number counts ex- 
hibit these tails at high density. Moreover, we verified that, de- 
populating randomly the Omocks to have them matching the 
observed number counts as a function of /-band, and then re- 
computing 6, does not suppress these tails. This means that in 
the Munich semi-analytical model, the 3D spatial distribution 
of galaxies is intrinsically different from that observed in the 
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Fig. 4. Density contrast distribution in the different redshift bins 
and luminosity limits explored. Red thick solid line: VVDS sam- 
ple; black dotted line: Omock with lowest number counts; black 
thin solid line: average of the Omocks; grey area: l-cr scatter of 
Omocks. 



real Universe, and that the excess of the clustering signal in the 
mocks is only in part due to an excess of low- to intermediate- 
mass galaxies. In Paper I, we suggested that this might be due 
to the assumption of a WMAP1 cosmology in the Millennium 
Simulation, and in particular to the use of a high normalisation 
of the power spectrum (cr 8 ). A lower value of cr 8 w ould reduce 
the o verall density contrast at any given redshift dWang et al.l 
12008b . Indeed, in Paper I we showed that, converting the cor- 
relation functions in the model to those expected assuming a 
lower value of <x 8 , the signal decreases at all scales. However, 
in Paper I we did not cha nge other cosmological parameters. 
Recently, iGuo et al" IJ20T2l) rescaled the Millennium simulation 
to the WMAP7 ( Komatsu et alJl20l"lh cosmological parameters 
values. They showed that the effects of the decreased erg are 
compensated by the higher value of O m . Therefore, the assump- 
tion of a WMAP1 cosmology cannot explain the excess of clus- 
tering in our mocks. 



4.2. The colour distribution 

C06 empirically defined red and blue galaxies to be the two 
extremes of the u*-g' colour distribution. Namely, red galaxies 
are defined as those with u* - g' > 1.1 and blue as those with 
u* - g' < 0.55. These colour cuts roughly correspond to the two 
peaks of the bimodal colour distribution in the VVDS. C06 kept 
these limits fixed for all luminosities and redshift bins consid- 
ered. 
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Fig. 5. Colour distribution in the different redshift bins and lumi- 
nosity limits explored. Red thick solid line: VVDS sample; black 
thin solid line: average of the Omocks; grey area: l-cr scatter of 
Omocks. Red short arrows represents the fixed colour cuts used 
to define red ((u* - g') > 1.1) and blue ((«* - g') < 0.55) galax- 
ies in C06. Blue long arrows indicates which colour cuts should 
be used in the mocks to have the same percentage of red and 
blue galaxies (irrespectively of environment) as in the VVDS, 
for each redshift bin and luminosity limit. 



From Paper I and from Sect. [3] we know that the colour 
distribution in mocks is different from the observed one. Fig. [5] 
shows the colour distributions in Omocks and VVDS for the dif- 
ferent redshift bins and luminosity limits considered in this pa- 
per. The l-cr scatter for the Cmocks is narrower than the Omocks 
one, but the average values of the two kinds of mocks are very 
close. We note that at intermediate redshifts (0.6 < z < 1.2), the 
blue cloud is more populated in the mocks than in the VVDS, 
especially for bright galaxies (Mb < -20). At these redshifts 
and luminosities, the blue peak of the bimodal colour distribu- 
tion in the mocks is from 30% to 100% higher then the VVDS 
one, while having roughly the same width. Still, if we consider 
only the tail of bluest galaxies (u* - g' < 0.55, as defined in 
C06), these are less in the mocks than in the VVDS, especially 
for faint galaxies. In contrast, the peak of the red galaxy pop- 
ulation is higher in the mocks than in data for fainter galaxies, 
and it corresponds to redder colours. These results are in agree- 
ment with those presented in Paper I, and with the discrepancies 
discussed in Sect. [3] 

Therefore the mix of galaxy populations (and colours) is dif- 
ferent in mocks and VVDS. As a consequence, if we used the 
same colour cuts in the mocks as in the VVDS, we will have 
different fractions of red and blue galaxies (independently of en- 
vironment). Fig.|6]shows the fraction of red and blue galaxies as 
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Fig. 6. Fraction of red ((«* - g') > 1.1, left panel) and 
blue ((m* - g') < 0.55, right panel) galaxies as a function 
of redshift, for different luminosity threshold. Lines are for 
VVDS data, triangles for Omocks and squares for Cmocks. 
The colour/thickness code is the following: from magenta/thin 
to blue/thick we consider galaxies with (Mb - 5 log h) < 
-19.0, -19.5, -20.0, -20.5, -21.0. Each triangle(/square) repre- 
sents the mean value of all the Omocks(/Cmocks). 



defined in C06 in the VVDS, as a function of redshift and for the 
different luminosity limits considered. The corresponding frac- 
tions in the mocks are different from those in the VVDS sample. 
Also, the colour dependency on luminosity is less evident in the 
mocks than in the VVDS. 

Fig. [6]also shows that the fraction of red galaxies is slightly 
lower in Omocks than in Cmocks at all redshifts, especially at 
z < 0.6, where it is lower by 5-8% depending on the luminos- 
ity cut. It is known that the VVDS observational strategy misses 
a very small fractio n (~ 4%) of very red galaxies in particu- 
lar at low z (see e.g. iFranzetti et al.ll2007b . This bias is mainly 
due to the optimisation for slit positioning, that targets prefer- 
entially smaller galaxies in /-band. At low redshift, brighter and 
bigger galaxies in the observed /-band are the early types. We 
will show that this small loss of red galaxies at z < 0.6 does not 
alter the colour-density relation in Omocks, and we believe that 
this applies to the observed data as well. We also verified that the 
global sampling rate does not vary as a function of the density 
computed on a 5/i~'Mpc scale in the Omocks, for the redshift 
bins and luminosity cuts considered. 

Fig. |5] also shows the two colour cuts of C06 (red arrows), 
and the colour cuts we should use in the mocks in order to have 
the same fractions of red and blue galaxies as in the VVDS 
for each redshift bin and luminosity cut (blue arrows). For this 
plot we are considering all galaxies, independently of their en- 
vironment. To compute the colour-density relation in mocks, we 
have taken the colour cuts corresponding to the same fraction 
of red and blue galaxies as in the VVDS, as this choice follows 
the same rationale used in C06, i.e. considering the extremes of 
the colour distribution. This makes the comparison with VVDS 
more straightforward, and provides the same overall normalisa- 
tion in the colour-density relation as in the VVDS. As a test, 
we have also computed the colour-density relation in the mocks 
using colour cuts that correspond roughly to the two peaks of 
the mock colour distribution. As for the VVDS, we kept these 
cuts constant at all redshifts and for all luminosities. The colour- 
density relations obtained in this way are very similar to those 
obtained using variable cuts that reproduce the observed frac- 
tion of red and blue galaxies. So our results are robust against 
this choice. 



4.3. The colour-density relation up to z ~ 1.5 

We reproduced the analysis on the colour-density relation as in 
C06, using both Cmocks and Omocks catalogues. 

Fig. [7] shows the fraction of red (f re d) an d blue (fblue) galax- 
ies as a function of the density contrast 5, for different redshift 
bins and for different luminosity thresholds. Red and blue sym- 
bols are for Omocks, while green and orange shaded areas show 
the error contours for the VVDS. Fig. |8]shows the corresponding 
quantities for Cmocks. The results from Omocks have larger er- 
ror bars due to larger Poisson error, but the average is very close 
to that of Cmocks. In each mock, we computed f rec i and fbi ue 
as a function of 6 in three equipopulated density bins (choosing 
the median 8 value as representative for the bin). The points in 
Fig. |7] and [8] are the average of the mocks, on both x- and y-axis. 
Vertical error bars are the sum in quadrature of the rms around 
the mean of the 23 mocks and the typical error (Poisson) of one 
single mock. In this way, we simultaneously account for cosmic 
variance among different realizations, and for the typical Poisson 
noise in a single realization. 

4.3.1 . The effects of a <1 00% sampling rate on the 
colour-density relation 

We find that the average colour-density relation in Omocks is 
very similar to the one found in Cmocks. This means that, on 
average, the VVDS observational strategy does not significantly 
alter the observed environmental trends for galaxy colours, on 
a ~ 5h _I Mpc scale. This does not exclude the possibility that 
there may exist at least one mock in which the colour-density 
relation is strongly affected by the VVDS observational strategy. 
In order to address this question, we performed a linear fit of f re d 
as a function of 5 in each luminosity and redshift bin considered, 
and compared the slope of this linear fit in Cmocks and Omocks. 
We did the same calculation for fu ue . 

In Fig.|9j we show the slopes obtained for each mock, in the 
four redshift bins considered and for galaxies with Mb < -20. 
(i.e., the third row of Fig. |7}- We plot the slope values for both 
Cmocks and Omocks. The x-axis value is arbitrary, but for each 
couple of Cmocks and Omocks extracted from the same light 
cone we use the same x-axis value for f rei j and fblue- Fig.[9]shows 
that a positive slope in Cmocks never becomes negative in the 
Omocks, or vice-versa, but for a couple of exceptions at 0.25 < 
z < 0.6 which may be due to low number statistics. The figure 
also show the average slopes for Cmocks and Omocks, together 
with the slopes from the VVDS sample. 

This Figure shows that going from Cmocks to Omocks, the 
colour-density relation becomes slightly shallower, but does not 
disappear nor reverse. This suggests that the trends found in C06 
(the flattening and possibly the reversal of the colour-density re- 
lation going to high redshift) are not due to biases in the VVDS 
observational strategy (uneven and < 100% sampling rate, field 
shape, etc.). 

4.3.2. Observed data versus mocks 

In the previous Section, we showed that the VVDS strategy does 
not significantly alter the colour-density relation potentially ob- 
servable in a survey with 100% sampling rate. Now we analyse 
the colour-density relation in the VVDS observed sample and the 
one found in mocks, with the aim to contrast galaxy evolution in 
simulations and in the real Universe. 

Some interesting trends are visible in Fig. Q and [8] First, 
the density contrast in mocks spans a larger range than ob- 
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Fig. 7. Colour-density relation in Omocks: blue squares represent 
the fraction of blue galaxies (fbiue) an d red triangles the fraction 
of red ones {freed- x- and y-axis positions are mean values of 
the x- and y-axis values for the 23 mocks. Vertical error bars in- 
cludes the scatter among the 23 cones and the typical (mean) 
error in each mock on the computation of fbiue or f re d- Orange 
and green contours represent f m t and f\ ue in the VVDS sam- 
ple (from C06). The numbers quoted in each panel are the total 
number of red and blue galaxies in the VVDS data and in the 
mocks (average on the 23 mocks), with the same colour code as 
for the symbols. For VVDS observed data we define red galaxies 
those with (u* - g') > 1.1 and blue those with (m* - g') < 0.55, 
at all z and for all luminosities. In the mocks, the definition of 
red and blue galaxies varies with z and luminosity, in order to 
have the same fraction of red and blue galaxies as in the VVDS 
(irrespectively of environment) in each redshift bin and for each 
luminosity limit (see the blue arrows in Fig. [5). 



served data, extending towards higher densities. This mirrors the 
stronger clustering found in mocks and discussed in Paper I and 
in Sec. 12.21 Second, error bars for Cmocks are smaller than those 
of VVDS, although they include the cosmic variance among the 
23 mocks. This is because the Poisson noise is much lower, as 
can be seen comparing the number of galaxies shown in the la- 
bels of Fig. |8] In contrast, error bars for Omocks are larger than 
those of VVDS, because their Poisson noise is similar but they 
include also cosmic variance. These trends are similar for all 
redshift bins and luminosity thresholds considered. 

At a fixed luminosity threshold, the colour-density relation in 
VVDS is steeper than in mocks at 0.25 < z < 0.6, very similar 
to that in mocks at 0.6 < z < 0.9, and shallower than in mocks 
at z > 0.9. In particular, the colour-density relation in the VVDS 
data is almost flat at 0.9 < z < 1 .2, and it seems to be inverted 
at z > 1.2 (i.e., at these high z blue galaxies reside preferentially 
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Fig. 8. Like in Fig. [7] but for Cmocks. 
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Fig. 9. Slopes of the colour-density relation linear fits, for the 23 
Cmocks and the 23 Omocks, for galaxies with Mb < -20, in the 
four redshift bins (label on the top), like in the third row of Fig. [7] 
and Fig. [8] Red and blue small diamonds: red and blue galaxies 
in each Omock; light purple and cyan small diamonds: red and 
blue galaxies in each Cmock. Red and blue big triangles: mean 
values of red and blue small diamonds, with the vertical error bar 
being the rms among all the 23 mocks. Light purple and cyan big 
squares: mean values of light purple and cyan small diamonds, 
with their rms. Orange and green circles: slope, and its error, 
for VVDS observed data (from C06), for red and blue galaxies 
respectively. 



in high density regions). In contrast, the mock colour-density 
relation varies only weakly as a function of redshift with no sig- 
nificant flattening, and definitely no inversion at higher redshift. 
We analysed in details the colour-density relation in the cou- 
ple of Cmock/Omock with the lowest number counts. As for the 
density distribution shown in Fig. [4] this mock shows a colour- 
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density relation very close to the average of the 23 mocks. We 
will discuss the implications of these general trends in Sect. [6] 

Fig. [10] shows the slopes obtained from fitting the colour- 
density relation in the mocks and compares them with those 
obtained from the VVDS data. Although the error bars for the 
VVDS data are large, the evolution of the colour-density rela- 
tion with redshift is clear. In contrast, no evolution is found in 
Omocks, nor in Cmocks. Moreover, while in the VVDS data 
there is a trend for steeper colour-density relation for brighter 
galaxies, this trend is much less clear (if any) in mocks. We note 
that the dependence of the colour-density relation on galaxy lu- 
minosity is a controvers ial issue in the liter ature. It has been 
found in C06, but not in ICooper et al.1 (12007 ), who studied th e 
colour-density relation in the DEEP2 survey ( Davis et aTll2003l) . 
selected with a sligh tly brighter flux limit than the VVDS. In 
ICucciati et al.l (1201 Oh . we did not fin d any luminosity depen- 
dency in the zCOSMOS Bright survey dLillv et alj|2009t) . which 
is much brighter {Iab ^ 22.5) than the sample used in C06. 

The comparison between the slopes measured in the mocks 
and those measured from data may suffer from the fact that the 
density distribution in mocks extends to higher densities than ob- 
served. Because of this, in almost each panel of Fig. [7] and [8] the 
point at highest density has a larger x-value in mocks than in the 
VVDS, even if f re< i or fbi ue are very similar. Thus the slope may 
be flatter in mocks than in the VVDS, only because of the differ- 
ent density distribution. To do a more fair comparison, we com- 
puted the fractional increment of f re d and fuue from the lowest 
to the highest density bin in each panel of Fig.|7]and[8] The re- 
sults, for both observations and mocks, are very similar to those 
of Fig. |T0j and confirm that there is no significant flattening of 
the colour-density relation in the mocks. 

The results discussed above confirm that, at least in the 
semi-analytical model used in our study, the environment shapes 
galaxy evolution (as a function of z) in a different way than we 
observe in the real Universe. 

5. Why is there a colour-density relation? 

Different methods to parameterise the local density around 
galaxies have been used in the literature for different galaxy sam- 
ples. Often it is the survey strategy itself that dictates the opti- 
mal (less biased) environment parameterisation for each specific 
data set. It is not yet clear which is the 'physical meaning' of 
these parameterisations, e.g. how the density field compares to 
the underlying dark matter distribution, and in particular how the 
estimated 'density' relates to the mass of the DM halos in which 
galaxies reside. 

Some recent studies (H aas et al 1 1201 2tlMuldrew et al.|[2012b 
have focused on a comparison between different environmen- 
tal definition and the information on the parent DM halo mass. 
These studies have considered results at z ~ 0, and have not 
entered the details of different observational strategies and/or se- 
lections. Given the evolution of the mean density and of the mass 
growth of structures, it would be very interesting to extend this 
detailed analysis to higher redshift. It should be noted, however, 
that results might well depend on the details of each particular 
survey and, as such, they are difficult to generalise. 

We have taken advantage of the available MILLENNIUM 
mocks to understand the origin of the colour-density relation 
as observed in VVDS data. To do this, we explored the re- 
lationship existing between the estimated density contrast and 
the mass of the DM haloes in which galaxies reside. For each 
galaxy in our mocks, we retrieved its parent DM halo us- 
ing the public database built for the Millennium Simulation 



(iLemson & Virgo Consortiumll2006l) . Here, a 'DM halo' corre- 
sponds to a halo identified in the Millennium Simulation using 
a friends-of-friends algorithm, with a linking length of 0.2 in 
units of the mean interparticle separation. Our results are shown 
in Fig. Q~TJ for galaxies with Mb < -20. We repeated our anal- 
ysis for both Cmocks and Omocks (top and bottom panels of 
the figure), in order to test the influence of the VVDS obser- 
vational strategy. In each panel, we distinguish central galaxies 
from satellite galaxies. On the right vertical axes of Fig. [TTJ we 
indicate the typical radius (in /t^'Mpc) of haloes with mass given 
by the corresponding value on the left axifl This shows that the 
scale on which 5 is computed is much larger than the size of the 
structures in which galaxies reside. 

Fig. [TTJ shows a very general trend: galaxies belonging to 
massive halos (mass > 10 b Mq//i) reside only in over-dense 

regions on a 5h~ l Mpc scale, even if the virial radius of such 
halos is much smaller than the filtering scale used. This hap- 
pens because the density in these regions is boosted to high val- 
ues by the large number of galaxies residing in these massive 
haloes, and is not affected significantly by the large scale struc- 
ture around them. In contrast, lower mass halos span the entire 
density range, because their density within the virial radius is 
not very high, so the density on larger scales depends also on the 
surrounding structures. 

From Fig. [TTJ it is clear why we should expect a colour- 
density relation on a 5h~ l Mpc scale. There is an increase of 
the parent halo mass with increasing 6 for satellite galaxies, al- 
though with a large scatter. This relation is driven by the large 
number of satellites in the most massive halos. The Figure shows 
that, if a satellite has a high measured 5 on the scale considered, 
it more probably belongs to a massive cluster than to a low mass 
halo. If we assume that clust ers have higher fractio n of red satel- 
lites than groups (see e.g. iDe Lucia et al.l [20T2h . we expect a 
colour-density relation for satellite galaxies. 

For central galaxies, the correlation between halo mass and 8 
on 5h~ l Mpc is weak, but there is still a trend for central galaxies 
in highest density regions being scattered to higher mass haloes. 
This happens because the most massive haloes tend to be clus- 
tered. Then, if we assume that central galaxies in more massive 
haloes are redder than central galaxies of lower mass haloes, we 
would expect to have a higher fraction of red galaxies in higher 
density regions. 

In Fig.[7]and[8] we showed that we do find a colour-density 
relation for model galaxies. We verified that this relation is main- 
tained when considering only centrals or only satellites galaxies. 
The colour-density relations for the two populations are similar, 
and in agreement with the global colour-density relation. 

As redshift increases, the VVDS flux limits select brighter 
and brighter galaxies. Thus, the fraction of observed satellites 
(that dominate the intermediate to faint end of the LF) decreases, 
so that the correlation between density and halo mass becomes 
less significant also for satellites. This is also evident comparing 
Cmocks with Omocks: the VVDS observational strategy applied 
in Omocks reduces the fraction of observed satellites in massive 
clusters, flattening the relation between halo mass and 5. This 
happens because the size of the slits in the spectrograph prevents 
us from targeting, in one single pass, galaxies with small pro- 
jected distances. The VVDS multi-pass strategy alleviates this 



4 The virial mass is computed using the simulated particles, as the 
mass enclosed within a sphere that corresponds to an overdensity of 
200 times the critical density. The virial radius is computed from the 
virial mass, through scaling laws based on simulation results and the 
virial theorem. 
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problem but small-scale very dense regions (like the central re- 
gions of galaxy clusters) are still undersampled with respect to 
regions that are less crowded. At the highest redshift considered, 
where the effect of the low sampling rate sums up with the effect 
of the flux limit, we do not see any clear correlation, even for 
satellite galaxies. 

In a survey like the VVDS, given its flux limit, the fraction 
of satellites galaxies (f sat ) is small, and it decreases with z. In 
particular, for galaxies with Mb < -20, it varies from ~ 15% to 
~ 8% going from the lowest to the highest redshift explored. So, 
we expect that central galaxies have an important role in shaping 
the observed colour-density relation. We will come back to this 
in Sec. [6] 



6. Discussion 

As shown in previous studies and above, the Munich semi- 
analytical model does not reproduce some of the galaxy proper- 
ties measured from the VVDS. In particular, in the model there 
is an excess of sources for Iab < 24, mainly due to an excess 
of faint red galaxies. The model also exhibits a slight excess of 
bright blue galaxies, and the galaxy clustering is overpredicted. 
In this paper, we have studied how in this model the 3D galaxy 
distribution affects galaxy colours up to z — 1.5. The aim of our 
analysis is to understand how environment affects galaxy evolu- 
tion modulating internal physical processes, and to what extent 
this is reproduced by the model used in this study. 

Our analysis demonstrates that the colour-density relation, 
computed on a 5/i~'Mpc scale, does not evolve significantly 
with redshift in the Munich semi-analytical model, at least from 
z = 0.25 to z = 1.5. This relation does not evolve even if we 
consider centrals and satellites separately. In contrast, significant 
evolution has been o bserved in different s amples and on differ - 
ent scales (e.g. C06, ICooper etal1l2007l ICucciati et alll2010h . 
and C06 even observed a possible inversion in such relation at 
z> 1.2. 

Given the striking difference between observational results 
and model predictions about the evolution of the colour-density, 
some questions arise: within a cosmological framework, do we 
expect an epoch when the colour-density relation is inverted? On 
which scale should such an inversion be measured? And are the 
physical processes responsible for such a relation included in the 
model that we have considered in this study? 

Qualitatively, we expect a reversal of the colour-density rela- 
tion, at an early epoch and on the scale of galaxy groups, due to 
starburst events in gas-rich mergers. In this scenario, two young 
and gas-rich galaxies will evolve differently if one remains in 
isolation (passive evolution) and the other merges with another 
gas-rich galaxy (this will trigger a starburst episode whose inten- 
sity depends on the mass ratio and on the amount of gas avail- 
able). During the merger, and for some time after it has been 
completed, the fraction of bright star-forming galaxies should be 
higher in high densities than in low densities. This simple sce- 
nario is complicated by the fact that lower density environments 
will also contain young star forming galaxies, although these 
might be fainter than the VVDS flux limit. Another complica- 
tion is due to dust attenuation: a large fraction of the galaxies that 
are classified as red might be forming stars at some significant 
rate, and this fraction might evolve as a function of cosmic epoch 
and/or luminosity and environment. For example. iDe Lucia et al.l 
(I2012h find that the fraction of red star forming galaxies de- 
creases with increasing mass and with decreasing distance from 
the cluster center. Moreover, it has been observed that dust at- 
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Fig. 11. Total halo mass as a function of the local density contrast 
computed in this work, in different redshift bins (as in the labels), 
for galaxies with Mb < -20. The values of the density contrast 
are indicated as log{\ + 6) in the bottom x-axis and as S in the 
top x-axis. The halo is the FOF halo in which the galaxy reside 
(the halo of the group/cluster of which the galaxy is member). 
Four top panels: Cmocks. Four bottom panels: Omocks. Red cir- 
cles and contours refer to central galaxies, blue squares and con- 
tours to satellite galaxies. Each panel includes all the 23 mocks. 
Lines represent the isodensity contours of all the galaxies in the 
panel. Not to crowd the plot, single galaxies are plotted as small 
dots only outside the lowest density contour. Filled red circles 
(/blue squares) are the median values, in density bins, for central 
(/satellite) galaxies. The error bars represents the 16% and 84% 
of the halo mass distribution in each density bin. On the right 
vertical axis of each panel, we show the typical virial radius (in 
/j~'Mpc) corresponding to the halo mass on the left vertical axis. 
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tenua tion clearly evolves with redshift (see e.g. ICucciati et alJ 
l2012h . 

A similar scenario , based on starburs t events in wet mergers , 
has been proposed bv lElbaz et al 1 (I2007h . who find an inversion 
of the star formation - density relation comparing the GOODS 
fields at high redshift with SDSS data. In particular, they find the 
mean star formation rate to be larger in high densities at z ~ 1 , on 
a 1.5/i~'Mpc scales. Their galaxy selection is similar to the one 
adopted in C06 at the same redshift (Mb < -20). In their sample, 
these galaxies are often located in correspondence of local den- 
sity enhancement on the scale of clusters/groups (~ l/r'Mpc). 
lElbazetal.ld20Q7l) also measure t he SF-density relation in S AMs 
(th e model by Croton et al.l2006l applied on the Millennium Run 
bv lKitzbichler & Whitdl2007ir applying the same methods as in 
their observational sample, and they did not find any reversal at 
z ~ 1, but a mild reversal of the SF-density relation at z ~ 2. 
I Wang et all (120071) carried out a similar analysis focusing on the 
relation between the average D4000A break and the local den- 
sity computed on a ~ 2/j~ 1 Mpc scale. In agreement with what 
discussed above, they find that in the Munich semi-analytical 
model there is no significant weakening of the D4000A-density 
relation up to redshift ~ 3. 

It is interesting that in C06 we found a possible reversal of 
the colour-density relation on ~ 5/j~'Mpc scales, that is much 
larger than those of galaxy groups. A few previous studies have 
argued that large scale environmental trends a re the residual of 
the trends observed on much smaller scal es (iKauffmann et al.l 
2004t iBlanton et al.1 l2006t ICucciati et all 1201 Ol and see also 
Wilman et al 1 12010ft . We can not establish if this is the case also 



CMOCKS 



in C06, as the scale used to compute the density contrast was 
the smallest allowed by the VVDS sampling rate to guarantee a 
reliable reconstruction of the density field (see the tests in C06). 
In addition, the ~ 5/i~'Mpc scale should be compared to the 
typical scales of processes happening 'around' clusters, like for 
instance the infall of galaxies and groups onto larger structures. 
One may also expect that such scale varies with redshift, as struc- 
tures grow with time. 

It is interesting that in the Munich semi-analytical model we 
find the same (not-evolving) colour-density relation for satellites 
and central galaxies separately (see Sec. |5]l. We should consider 
that some physical processes that might be at play in determin- 
ing the evolution and the inversion of the colour-density relation 
are modelled in a rather crude way. For example, starbursts trig- 
gered by mergers are 'instantaneous' (i.e. their time-scale cor- 
responds to the integration time-scale adopted in the model), 
and the model does not account for star formation episodes that 
might be triggered during fly-by. In addition, in this model (like 
in most of the recently published ones) there is a significant ex- 
cess of red (passive) galaxies with respect to observations. In 
particular, the model galaxies in the region of infall around clus- 
ters are redder than in the real Universe, because most of the 
galaxies in a group that is infalling in a cluster will be already 
too red in the model. 

Since the model does not match the expected evolution of the 
colour-density relation as a function of redshift, we can modify 
the colour distribution in the model by making simple assump- 
tions, and see how the colour-density relation changes. For ex- 
ample, we verified what happens if we assume that the colour 
of central galaxies is correctly reproduced by the model, and 
that satellites are either all red or all blue. The results of this 
test are shown in the top panel of Fig. [12] Since f sat increases 
with density, at all z, if satellites are all red the increase of f re d 
as a function of density would be even steeper than what found 
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Fig. 12. Top panel. Slope of the linear fit of the /, e j-density 
relation (thick red lines) and of the ft,i„ e -density relation (thin 
blue lines) as a function of redshift in Cmocks. Only galaxies 
withMfi < -20 are considered. Solid line: slopes as found in 
the model, the same as in Fig. [10] Dashed and dot-dashed lines: 
slopes for two different assumptions about the colour distribu- 
tions of satellite galaxies. We assume that centrals galaxies have 
the same colour as in the model, but now we assume that satel- 
lites are either all red (dashed line) or all blue (dot-dashed line). 
Bottom panel. Slope of the linear fit of the / rK /-density relation 
as a function of redshift in Cmocks, for three different assump- 
tions about the colour distributions of central and satellite galax- 
ies. Solid line: all satellites are red, and all centrals are blue. 
Dotted line: all satellites and centrals residing in DM halos with 
virial mass > 10 13 M G //z are red, while all other centrals are blue. 
Dashed line: like the dotted line, but in this case the red centrals 
are those residing in DM halos with virial mass > 10 12 M o //j. 
For this plot we use the 23 Cmocks together, and we keep f sat 
increasing with 6 and decreasing with z as we find in the mocks. 



in the model. So this would increase the disagreement between 
data and models. If we assume that satellites are all blue, which 
is unrealistic, we find that the f rec [-5 relation flattens at all red- 
shifts, with the flattening being a bit more significant at low red- 
shift. This is because f sal decreases with increasing redshift (see 
Fig. [Til . This is, however, not enough to have a flat relation at 
z > 1.2, as observed for VVDS. In addition, assuming all satel- 
lites are blue, would give an almost flat fbiue-5 relation at z > 1.2 
but would invert the relation at z < 0.9, contrary to what is ob- 
served. 

This suggests that an inaccurate modelling of the physics of 
satellite galaxies cannot be solely responsible for the disagree- 
ment we find between models and data, and that such a dis- 
crepancy also signals the need for a revised treatment of central 
galaxies. 

In order to gain insight on this, let us assume that the model 
predicts the correct increase (decrease) of f sat as a function of 
density (redshift), and that all satellites are red while all centrals 
are blue. Since f sar decreases with increasing redshift, this would 
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make the f re d-S relation flatten with increasing redshift, as shown 
in the bottom panel of Fig. [12] The figure shows that a flattening 
would be observed also in the case we assume all satellites and 
central galaxies of haloes more massive than 1O 13 M //1 (dotted 
line) are red, and that all other centrals are blue. On the contrary, 
if we assume to be red all satellites and all central galaxies of 
haloes more massive than 10 l2 M Q /h (dashed line), the evolution 
of the slope of the f re d-5 is not monotonic. 

The toy-model considered is simplistic, because for instance 
we consider all satellites to be red, and we do not include 
any redshift evolution of the colour-mass relation for centrals. 
Anyway, this simple exercise tells us that, for a survey like 
VVDS, the role of central galaxies in the evolution of the colour- 
density relation is very important, and that also the modelling of 
these galaxies needs to be revised in the model used in this study. 

Further work is therefore needed in order to understand what 
are the physical processes driving the flattening (and possibly 
inversion) of the colour-density relation in real observations, 
and how to include them properly in the context of hierarchi- 
cal galaxy formation models. An interesting approach would be 
to use the accretion h istories of model galax ies in different den- 
sity bins, as done in iDe Lucia et al.l d2012l) . If galaxies in low 
and high density regions have spent different fractions of their 
life-time in high density environments, we expect them to have 
different colours. Using this approach, it will be possible to com- 
pare, and possibly to link, the formation histories of galaxies at 
low and high redshift, and to shed light on the evolution of the 
colour-density relation. We will address this issue in future work. 

7. Summary and Conclusions 

In this paper we used galaxy mock catalogues constructed from a 
semi-analytic model coupled to a large cosmological simulation, 
to compare t he observed colour-d ensity analysis up to z ~ 1.5 
presented in ICucciati et al.l (120061) with model predictions, re- 
producing carefully in the mocks the observational selection and 
strategy adopted for the VVDS Deep sample. 

In particular, we used 23 galaxy mock catalogues, with the 
same flux limits adopted in the VVDS-Deep survey (Cmocks). 
The mocks are extracted from the Mill ennium Simulation, 
with a pplied the semi-analytical model by IDe Lucia & Blaizotl 
(I2007h . From each of these mocks, a sub-sample of galaxies has 
been extracted (Omocks) mimicking the VVDS survey obser- 
vational strategy (sampling rate, slit positioning, etc.). We then 
computed the galaxy luminosity function and the colour-density 
relation in the mocks, using the same methods employed for the 
VVDS survey. Our results can be briefly summarised as follows: 

1) the mocks based on the Munich semi-analytical model con- 
tain on average ~ 10% more galaxies than the VVDS sam- 
ple. This is consistent with results obtained in Paper I. 

2) The rest-frame B-band LF in mocks is in rough agreement 
with the one derived from the observations, but it has a 
slightly steeper slope, especially at 0.2 < z < 0.8. This is 
expected, given the excess of galaxies in the mocks in this 
redshift range. We have also computed the LF for early and 
late type galaxies separately, and shown that the model over- 
predict the number densities of faint early type galaxies and 
that of bright late type galaxies. 

3) The density distribution computed with the same method as 
in the VVDS has more prominent tails towards higher den- 
sities in the mocks, at any redshift and luminosity explored. 
This is not due to the larger number counts in mocks, but 
it is related to an intrinsically different galaxy distribution, 



which also reflects in a stronger clustering signal in mocks 
(see Paper I). 

4) The colour-density relation in Omocks is in very good agree- 
ment with the one in Cmocks (being only more noisy and 
slightly less significant in Omocks). This enhances the con- 
fidence that the evolutionary trend of the colour-density re- 
lation in the VVDS-Deep survey is not caused by any severe 
bias introduced by the survey observational strategy. 

5) The colour-density relation in mocks does not evolve signifi- 
cantly from z = 0.25 to z = 1.5, in contrast with a significant 
evolution measured in the data. In particular, we find no flat- 
tening (or inversion) of the colour-density relation at higher 
redshift, at least up to the redshift explored by the VVDS 
data and on the same scale. 

6) Given the relation between the measured density contrast 
and the virial mass of the halos where galaxies reside, we do 
expect and find a colour-density relation both for central and 
satellites galaxies. Both are very similar to the relation found 
for the global population. We argue that the lack of evolution 
in the colour-density relation in mocks cannot be due only to 
inaccurate prescriptions for the evolution of satellites galax- 
ies, and that also the treatment of the central galaxies has to 
be revised. 

A reversal of the colour-density relation is expected in a sce- 
nario in which wet mergers of young galaxies trigger an en- 
hancement of star formation in the interacting galaxies. In this 
scenario, a reversal of the colour-density relation should be ob- 
served on the scale of galaxy groups. A reversal of the star 
format ion-density relatio n has been observed on such scales at 
z ~ 1 (lElbaz et al.l 120071) . However, it is not straightforward to 
correlate galaxy SF with colours, because of dust reddening. So 
it is not clear that what we measure on a 5/i~'Mpc scale is a 
mirror of what is happening on smaller scales. The lack of evo- 
lution of the colour-density relation in the model, on a 5/z~'Mpc 
and at least up to z = 1.5, suggests that the evolution has hap- 
pened at z > 1.5, and/or that model galaxy colour is affected 
by environment on scales much smaller than 5/i~'Mpc, with no 
corresponding trends on larger scales. 

The disagreement between the evolution of the colour- 
density relation in mocks and in observed data deserves further 
investigation as it can clarify what are the physical processes 
driving the observed flattening (and inversion), and to what ex- 
tent these physical processes are included in recent models of 
galaxy formation. It would be also important to disentangle the 
role of central and satellites galaxies. In this work we focused 
on mocks reproducing the VVDS observational strategy. With 
these mocks it is not possible to go further in this analysis, be- 
cause the flux limits prevent us from having enough galaxies to 
study in details the small-scale environment at the highest red- 
shift explored, even if we had a 100% sampling rate. With mock 
catalogues it is possible to go beyond the luminosity and red- 
shift range explored by the VVDS, but currently there would 
not be a suitable counter-part among the available data sets. For 
example, iKnobel et alJ (1201 2l) successfully separate central and 
satellites galaxies in the zCOSMOS group catalogue, but this 
catalogue reaches only z = 1, and the flux limit {I a b < 2 2.5) 
is brighter than the one in the VVDS. iGerke et all (1201 2l) use 
data from the DEEP2 survey to compute a group catalogue that 
extends at z > 1, but they do not distinguish between central 
and satellite galaxies. Moreover, the DEEP2 flux limit is very 
similar to the VVDS one, so at z > 1 the study of satellite galax- 
ies would be difficult. Therefore, deeper spectroscopic surveys, 
with higher sampling rate, would be needed in order to shed light 
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on the physical processes establishing the observed relation be- 
tween colour and density. 

Acknowledgements. We would like to thank the referee for helpful suggestions 
that improved the paper. The Millennium Simulation databases used in this paper 
and the web application providing online access to them were constructed as part 
of the activities of the German Astrophysical Virtual Observatory. We are grate- 
ful to Gerard Lemson for setting up an internal data base that greatly facilitated 
the exchange of data and information needed to cany out this project, and for 
his continuous help with the database. OC thanks the INAF-Fellowship program 
for support. GDL acknowledges financial support from the European Research 
Council under the European Community's Seventh Framework Programme 
(FP7/2007-2013)/ERC grant agreement n. 202781. Part of this work was sup- 
ported by PRIN INAF 2010 'From the dawn of galaxy formation to the peak of 
the mass assembly'. The results presented in C06 have been obtained within the 
framework of the VVDS consortium, and we thank the entire VVDS team for its 
work. 



References 

Abbas, U. & Sheth, R. K. 2005, MNRAS, 364, 1327 

Bielby, R., Hudelot, P., McCracken, H. J., et al. 201 1 , ArXiv: 1111 .6997 

Blaizot, J., Wadadekar, Y, Guiderdoni, B„ et al. 2005, MNRAS, 360, 159 

Blanton, M. R., Eisenstein, D., Hogg, D. W„ & Zehavi, I. 2006, ApJ, 645, 977 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Brazual A., G. & Chariot, S. 1993, ApJ, 405, 538 

Chabrier, G. 2003, PASP, 115, 763 

Coleman, G. D., Wu, C, & Weedman, D. W. 1980, ApJS, 43, 393 

Cooper, M. C, Newman, J. A., Coil, A. L., et al. 2007, MNRAS, 376, 1445 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 1 12, 839 

Croton, D. J., Springel, V, White, S. D. M., et al. 2006, MNRAS, 365, 11 

Cucciati, O., Iovino, A., Kovac, K., et al. 2010, A&A, 524, A2 

Cucciati, O., Iovino, A., Marinoni, C, et al. 2006, A&A, 458, 39 

Cucciati, O., Tresse, L., Ilbert, O., et al. 2012, A&A, 539, A31 

Davis, M., Faber, S. M., Newman, J., et al. 2003, in Discoveries and Research 

Prospects from 6- to 1 0-Meter-Class Telescopes II. Edited by Guhathakurta, 

Puragra. Proceedings of the SPIE, Volume 4834, pp. 161-172 (2003)., ed. 

P. Guhathakurta, 161-172 
de la Torre, S., Meneux, B., De Lucia, G., et al. 2011, A&A, 525, A125 (Paper 

I) 

De Lucia, G. & Blaizot, J. 2007, MNRAS, 375, 2 

De Lucia, G., Kauffmann, G., & White, S. D. M. 2004, MNRAS, 349, 1 101 
De Lucia, G., Springel, V, White, S. D. M., Croton, D., & Kauffmann, G. 2006, 
MNRAS, 366, 499 

De Lucia, G., Weinmann, S., Poggianti, B. M., Aragon-Salamanca, A., & 

Zaritsky, D. 2012, MNRAS, 423, 1277 
Efstathiou, G., Ellis, R. S., & Peterson, B. A. 1988, MNRAS, 232, 431 
Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 33 
Farouki, R. & Shapiro, S. L. 1981, ApJ, 243, 32 

Fontanot, F., De Lucia, G., Monaco, P., Somerville, R. S., & Santini, P. 2009, 

MNRAS, 397, 1776 
Franzetti, P., Scodeggio, M., Garilli, B., et al. 2007, A&A, 465, 711 
Gavazzi, G., Pierini, D., & Boselli, A. 1996, A&A, 312, 397 
Gerke, B. F, Newman, J. A., Davis, M., et al. 2012, ApJ, 751, 50 
Gunn, J. E. & Gott, J. R. I. 1972, ApJ, 176, 1 
Guo, Q., White, S., Angulo, R. E., et al. 2012, ArXiv: 1206.0052 
Haas, M. R., Schaye, J., & Jeeson-Daniel, A. 2012, MNRAS, 419, 2133 
Hatton, S., Devriendt, J. E. G., Ninin, S.. et al. 2003, MNRAS, 343, 75 
Ilbert, O., Tresse, L., Arnouts, S., et al. 2004, MNRAS, 351, 541 
Ilbert, O., Tresse, L., Zucca, E., et al. 2005, A&A, 439, 863 
Iovino, A., McCracken, H. J., Garilli, B., et al. 2005, A&A, 442, 423 
Kaiser, N. 1987, MNRAS, 227, 1 

Kauffmann, G. & Haehnelt, M. 2000, MNRAS, 311, 576 

Kauffmann, G., White, S. D. M., Heckman, T. M., et al. 2004, MNRAS, 353, 
713 

Kitzbichler, M. G. & White, S. D. M. 2007, MNRAS, 376, 2 
Knobel, C, Lilly, S. J., Iovino, A., et al. 2012, ApJ, 753, 121 
Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 
Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 
Le Fevre, O., Mellier, Y, McCracken, H. J., et al. 2004, A&A, 417, 839 
Le Fevre, O., Vettolani, G., Garilli, B., et al. 2005, A&A, 439, 845 
Lemson, G. & Virgo Consortium, t. 2006, ArXiv:astro-ph/0608019 
Lilly, S. J., LeBran, V, Maier, C, et al. 2009, ApJS, 184, 218 
Lynden-Bell, D. 1971, MNRAS, 155, 95 

McCracken, H. J., Radovich, M., Berlin, E., et al. 2003, A&A, 410, 17 
Mo, H. J. & White, S. D. M. 1996, MNRAS, 282, 347 



Moore, B., Katz, N., Lake, G., Dressier, A., & Oemler, A. 1996, Nature, 379, 
613 

Muldrew, S. I., Croton, D. J., Skibba, R. A., et al. 2012, MNRAS, 419, 2670 
Peebles, P. J. E. 1980, The large-scale structure of the universe, ed. Peebles, 
P. J. E. 

Radovich, M., Arnaboldi, M„ Ripepi, V, et al. 2004, A&A, 417, 51 
Sandage, A., Tammann, G. A., & Yahil, A. 1979, ApJ, 232, 352 
Schechter, P. 1976, ApJ, 203, 297 
Schmidt, M. 1968, ApJ, 151, 393 

Springel, V, White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 
Springel, V, White, S. D. M., Tormen, G., & Kauffmann, G. 2001, MNRAS, 
328, 726 

Stringer, M. J.. Benson, A. J., Bundy, K, Ellis, R. S., & Quetin, E. L. 2009, 

MNRAS, 393, 1127 
Temporin, S., Iovino, A., Bolzonella, M„ et al. 2008, A&A, 482, 81 
Toomre, A. & Toomre, J. 1972, ApJ, 178, 623 

Wang, J., De Lucia, G, Kitzbichler, M. G., & White, S. D. M. 2008, MNRAS, 
384, 1301 

Wang, L., Li, C, Kauffmann, G., & De Lucia, G. 2007, MNRAS, 377, 1419 
Weinmann, S. M., Kauffmann, G., von der Linden, A., & De Lucia, G. 2010, 

MNRAS, 406, 2249 
Weinmann, S. M., van den Bosch, F. C, Yang, X., et al. 2006, MNRAS, 372, 

1161 

Wilman, D. J., Zibetti, S., & Budavari, T. 2010, MNRAS, 406, 1701 
Zucca, E., Ilbert, O., Bardelli, S., et al. 2006, A&A, 455, 879 



INAF-Osservatorio Astronomico di Trieste - Via Tiepolo 11, 1-34143 
Trieste, Italy 

2 INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 
Bologna, Italy 

3 INAF-Osservatorio Astronomico di Brera, Via Brera 28, 1-20021, 
Milan, Italy 

4 SUPA, Institute for Astronomy, University of Edinburgh, Royal 
Observatory, Blackford Hill, EH9 3HJ Edinburgh, UK 

5 Centre de Recherche Astrophysique de Lyon, UMR 5574, Universite 
Claude Bernard Lyon-Ecole Normale Superieure de Lyon-CNRS, 
69230 Saint-Genis Laval, France 

6 INAF-IASFBO, Via P. Gobetti 101, 1-40129, Bologna, Italy 

7 Aix Marseille Universite, CNRS, LAM (Laboratoire 
d' Astrophysique de Marseille) UMR 7326, 13388, Marseille, 
France 

8 The Andrzej Soltan Institute for Nuclear Studies, ul. Hoza 69, 00-681 
Warszawa, Poland 

9 Astronomical Observatory of the Jagiellonian University, ul. Orla 
171, PL-30-244, Krakow, Pola 



13 



O. Cucciati et al,; The VVDS colour-density relation at z < 1.5: model vs observations. 



Millennium 001b - UBVRI / all galaxies 
o^_L:r0.20-0.40 ' 



red=VVDS / black=001b 

r 



/ magenta=022b / cyan=036d 




-15 -20 
M Bi _-51og(h) 
Millennium 001b - UBVRI / CWW type 1 



21.0 -20.8 -20.6 -20.4 -21.0 -20.8 -20.6 -20.4 
M B-j -51og(h) 
Millennium 001b - UBVRI / CWW type 4 



•1 -i 



-5 - 



: 1 1 1 1 J 

r 0.20-0.40 - 


.1 1 1 | 1 1 1 | 1 1 1 | 1 1 1 | 1 1 1 | 1 1. 
r 0.40-0.60 J - 


■ r \\ ■ 

'^ L A — 1 — | — 1 — 1 — 1 — H-| — 1 — K 


- /H, i 

'p 1 v ^ ^ 

:■ In ■ 1 ' 1 ' 1 ■ 1 ■ h 1 ' | ' 1 ' 1: 


7O.6O-O.8O J -_ 
I 1 i 


^0.80-1.00 -_ 




: ( K : 


-20 -15 


-20 -15 



S 



•e- 

w 
o 




M B -51og(h) 



-15 

M. -51og(h) 



Fig. 3. Evolution of the luminosity function of galaxies in the mock samples, compared to the results from the VVDS. Upper left 
panel: Evolution of the luminosity function in the B-band for all galaxies in the Cmock 001b. Each panel refers to a different redshift 
bin, which is indicated in the label. The vertical dashed line represents the faint absolute limit considered in the S TY estimate. The 
luminosity functions are estimated with different methods (see text for details) but for clarity we plot only the results from C + 
(symbols) and STY (lines). Red filled circles and lines are used for observational measurements while black empty squares and 
lines are used for the corresponding measurements from the mock. Upper right panel: 68% confidence ellipses for the a and M* 
parameters obtained using all galaxies in mock 001b (black ellipse), mock 022b (magenta ellipse), and mock 036d (cyan ellipse). 
The red thicker ellipse corresponds to the VVDS measurements; the red thick lines in the first and the last redshift bins show the 
uncertainties on one of the parameters obtained keeping the other fixed. Lower left panel: Evolution of the luminosity function in 
the B-band for type 1 galaxies in the mock 001b. Lines and symbols have the same meaning as in the upper left panel. Lower right 
panel: Evolution of the luminosity function in the B-band for type 4 galaxies in the mock 001b. Lines and of symbols have the same 
meaning as in the upper left panel. 
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Fig. 10. Slopes of the linear fits of the colour-density relation, in each panel of Fig. |7]and Fig. [8] for f rec i on the left panel and fi,i ue 
in the right panel. Lines and shaded areas: VVDS observed data with error. Triangles: Omocks. Squares: Cmocks. Different colours 
of lines and symbols are for the different luminosity thresholds, as indicated in the labels. In each redshift bin, delimited by vertical 
dashed lines, all the points should be considered at the same (central) x-value, but they are shifted for clarity. 
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